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ABSTRACT 

The  purpose  of  this  thesis  'i!^^6- structure  a  periodic  review  model 
for  a  reparable  item  system  experiencing  normally  distributed  random 
demand.      The  items  are  issued  upon  demand,    repaired   (if  possible)    after 
use  or  failure,    and  subsequently  reissued.      The  model  addresses  the 
requirement  of  system  replenishment  as  items  are  lost  or  determined  to 
be  beyond  economical  repair.      The  inventory  of  ready-for-issue  items 
is  treated  as  two  separate  inventories,    i.  e.  ,    an  inventory  of  those  items 
received  directly  from  the  manufacturer  and  an  inventory  of  those  items 
received  from  the  repair  facility.     Annual  system  cost  expressions  are 
developed  as  a  function  of  the  desired  protection  level  and  length  of 
review  period. 
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A  -     Fixed  procurement  order  cost   (per  order)  . 

A  -      Fixed  repair  order  cost   (per  order)  . 

R 
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1.       INTRODUCTION 

An  existing  reparable  inventory  system  within  the  Naval  Aviation 
Supply  System  is  an  exceedingly  complex  probabilistic  system  involving 
tens  of  thousands  of  line  items  and  a  network  of  stock  points  and  repair 
points.     Anyone  attempting  to  model  a  system  of  this  complexity  is 
immediately  aware  that  the  model  can  never  be  more  than  a  partial 
representation  of  reality.      Despite  this  awareness  of  incompleteness, 
the  fact  that  models  have  proven  to  be  invaluable  aids  in  managing  com- 
plex systems  stands  as  a  primary  motivating  force  behind  the  efforts 
to  model  a  reparable  inventory  system. 

The  models  developed  by  Hoekstra   [1]  ,     McNall  and  Hatchett   [2]  , 
Schrady   [3]  ,    and  Allen  and  D'Esopo   [4]    attest  to  the  recent  efforts 
to  model  this  system.      The  purpose  of  this  paper  is  to  add  to  this  grow- 
ing knowledge  base  by  structuring  a  periodic  review  model  for  a  repa- 
rable item  system  experiencing  normally  distributed  random  demands. 
The  model,    in  turn,    provides  a  means  of  understanding  the  parameter 
interactions  and  cost -protection  trade-offs  that  exist  within  a  system  of 
this  nature.      The  periodic  model  techniques  used  by  Hanssman   [5]    in 
the  development  of  a  multilevel  production  control  model  are  relied  upon 
during  the  initial  phase  of  the  model  development. 

In  addition  to  the  model's  basic  assumptions  of  periodic  review  and 
normal  demand,    it  was  necessary  to  further  simplify  the  reparable  item 
system  by  making  additional  assumptions.      The  major  assumptions  are 
included  in  the  following  summary  of  the  scope  of  the  model. 


The  model  addresses  the  reparable  item  system  on  a  single  item 
basis.      The  system  is  considered  to  be  made  up  of  two  subsystems: 
(1)    the  repaired  item  subsystem,    and   (Z)    the  purchased   (new)    item 
subsystem.      The  purchased  item  subsystem  is  made  up  of  a  single 
manufacturer,    the  outstanding  purchase  order  quantities,    and  the  pur- 
chased item  portion  of  the  ready-for-issue    (RFI)    inventory  under  the 
control  of  the  inventory  control  point    (ICP)  .      The  repaired  item  sub- 
system is  made  up  of  a  single  repair  facility   (with  an  inventory  of  non- 
ready-for -issue    (NRFI)    items)  ,    outstanding  repair  order  quantities, 
and  the  repaired  item  portion  of  the  RFI  inventory  under  the  control  of 
the  same  ICP.      The  total  demand  each  period  is  apportioned  into  two 
segments  based  on  the  recovery  rate    r  ;    i.  e.  ,    r   percent  of  the  demand 
is  charged  against  the  ICP  repaired  item  inventory  and  the  remaining 
portion  is  satisfied  from  the  ICP  purchased  item  inventory.      The  as- 
sumption is  made  that  the  recoverable  portion  of  the  items  demanded 
each  period  will  be  returned  to  the  repair  facility  in  a  given  number  of 
periods  after  the  period  in  which  they  were  demanded. 

The  basic  model  assumes  that  repair  and  purchase  orders  are 
initiated  each  period  by  the  ICP.     A  sufficient  quantity  is  ordered  each 
period  to  bring  the  subsystems'  inventory  positions  up  to  established 
limits.      Separate  net  inventory  distributions  are  developed  for  the  ICP 
repaired  item  inventory,    ICP  purchase  order  inventory,    ICP  combined 
purchased  and  repaired  item  inventory,    and  the  repair  facility  inventory. 
The  relationship  between  the  high  limits  of  the  subsystem  and  the 
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corresponding  protection  level,    i.  e,  ,    one  minus  the  probability  of 
stockout,    is  determined  and  corresponding  system-cost  equations  are 
derived. 

The  equations  of  the  basic  model  are  then  slightly  modified  to  allow 
for  extension  of  the  repair  and  purchase  review  periods  in  multiples  of 
the  basic  review  period.      The  extension  does  not  require  that  these  two 
review  periods  be  of  equal  length,    as  was  assumed  in  the  basic  model. 
The  annual  operating  cost  trade-off  that  exists  between    (1)   period 
review  and  order  costs  and   (Z)    ICP  inventory  holding  costs  for  a  given 
protection  level  is  examined.     A  method  for  determining  the  operating 
procedures  that  minimize  the  expected  annual  operating  cost,    within  the 
evaluation  restrictions  of  the  extended  model,    are  then  investigated. 
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2.       A  PERIODIC  REPARABLE  MODEL 

2.  1    Description  of  Reparable  Item  System 

Figure   1    illustrates  the  periodic  structure  and  relative  leadtime 
relationships  of  the  reparable  item  system  addressed  and  modeled  by 
this  paper. 

The  description  of  the  system  will  commence  at  the  ICP.      The 
inventory,    I,    under  the  control  of  the  ICP  is  looked  upon  as  two 
inventories;    i.  e.  ,    those  items  received  directly  from  the  manufacturer, 
I     ,     and  those  items  received  from  the  repair  facility,    I     .      The  total 
withdrawal  quantity   (demand)    during  period   i   is  a  random  variable 
denoted  by   X.  .      It  is  assumed  that  period  demands,    X    ,    X    ,    .  .  .  ,    are 

independent  and  identically  distributed  as  a  normal  random  variable    X 

2 
having  mean  |j,       and  variance   o        .      The  recovery  rate,    r  ,    is  a  known 
X  X 

percentage;    that  is,     r   percent  of  the  items  demanded  each  period  will 
be  returned  in  a  reparable  state  to  the  repair  facility.     Each  period 

rX.    of  the  total  items  demanded  are  satisfied  from    I_    and  the  remain- 

X  2 

ing  portion,    (1    -  r)  X.   ,    of  the  total  demand  is  satisfied  from    I     . 

At  the  end  of  period  zero,    the  combined  ICP  inventory  position,    1(0)  , 
consists  of    I    (0)    plus    L(0)  .      The  outstanding  orders  at  this  moment 
in  time  are  shown  in  Figure  1  .      The  repair  order  quantity    q.    and  the 
purchase  order  quantity   P.    go  into  ICP  inventory  at  the  beginning  of 
period   i  .      The  purchase  order  leadtime  is    L   periods  and  the  repair 
order  leadtime  is    R   periods,    where  leadtime  is  defined  as  the  fixed 
interval  of  periods  between  the  time  the  ICP  decides  to  place  an  order 
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and  the  time  the  ordered  quantity  is  placed  in  the  ICP  controlled  inventory. 
Just  prior  to  the  beginning  of  period  1  ,    the  new  orders,     Q  and 

rv.  +    1 

P^        ,    ,    must  be  determined. 
L  +  1 

It  is  assumed  that  at  the  beginning  of  each  period  the  repair  facility 
receives  the  reparable  portion  of  the  system  demand  experienced  by 
the  system   K   periods  earlier,   where    K   is  defined  as  the  time  in 
periods  from  when  an  item  is  demanded  and  the  replaced  item  is  re- 
turned to  the  repair  facility.     An  ICP  repair  order    Q    is  received  by  the 
repair  facility  at  the  beginning  of  each  period.      Upon  receipt  of  the 
order,    the  repair  facility  immediately  commences  a  batch  repair  to 
fill  the  order  with  the  existing  stock  on  hand.      If  there  is  insufficient 
stock  on  hand  to  meet  the  order  at  this  point  in  time,    the  unfilled  portion 
of  the  new  order  is  included  in  the  batch  repair  initiated  at  the  beginning 
of  the  next  period. 

2.  2    ICP  Purchased  Inventory  Distribution  and  Order  Rules 

At  the  end  of  each  period,    a  new  purchase  order,     P  ,     is 

determined  by  the  ICP  and  transmitted  to  the  manufacturer.      It  is 
assumed  that  the  manufacturer  has  sufficient  raw  material  on  hand  to 
commence  production  of  the  ordered  items  immediately  upon  receipt  of 
the  ICP  order.      Therefore,    the  first    1      inventory  level  that  the  ICP  can 
control  is    I     (L  +   1  )  ,    where 

I^(L+1)    =     1^(0)    .    P^    +    P^    ,     ...     ,    P^    ,    p^^^ 


-    ((1    -  r)X^     +     ...     +     (1    -  r)  X^^  ^   )      .  (1) 
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Let   H      be  defined  as  the  high  limit  of  the  purchased  item  system,    which 
is  a  fixed  value  to  be  determined  by  management.     This  limit  represents 
the  maximum  level  of  items  in  the  purchased  item  system,    i.  e.  ,    ICP 
purchased  item  inventory  and  outstanding  purchase  order  quantities. 
Each  period  an  order  is  placed  to  bring  the  purchased  item  system 
back  up  to  this  level.     It  follows  that 

1=1 

or 


P  ,     =     new  order     =    H,      -     I    (0)     -LP. 

L  +   1  1  1  .    ,       1 

1=1 


Based  on  equation  (2)  ,    we  can  write  (1)    as  follows: 

L+1 
I    (L  +   1)    =    H       -      S     (1    -  r)X.      .  (3) 

i=l  ^ 

As  shown  in  Appendix  A,     I      is  a  normally  distributed  random  variable 

with  parameters 

p.^     =    H^     -    (L+   1)  (1   -  r)^^  (4) 

and 

o^    =    (L+  1)  (1   -  r)^a^  ,  (5) 

where  a  negative  inventory  is  considered  to  be  an  outstanding  backorder 
quantity  when  the  on-hand  inventory  level  is  zero. 

The  value  of  the  mean  parameter   p,       can  also  be  written  as  a  function 
of  a  desired  inventory  protection  level.      Given  a  desired  protection 
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level   a  ,    the  corresponding  probability  of  stockout  is     1    -  a   ,    which  can 
be  written 

P  [I      ^    0]     =    1  -  a     . 

This,    in  turn,    can  be  written  in  the  standard  normal  distribution  function 
form,    f  (i),    as  follows: 


0     - 

■      I     X 

1   -  a 
1 


U      -     M, 


or 


\     rr 

1 


Defining    N(a)    as  the  value  obtained  from  a  standard  normal 
distribution  table  corresponding  to  the  area   oi  under  the  curve, 

N(a)    =  -^       .  (6) 

^1 


Rearranging  and  replacing  a      with  its  equivalent  form  shown  in 


(5)  ,   this  can  be  written 


^.^     =    N(a)  (L+  1)^^^   (1   -  r)  o^       •  (7) 


From  equations  (7)  and  (4),    the  value  of  the  high  limit  of  the  pur- 
chase system   H      can  be  expressed  as  a  function  of  the  desired 


protection  level   a  and  the  parameters    p,        and  a^   • 

Hj     =    (1   -  r)   (N(a)  (L+  D^/^a^    +    {L+  l)^,^)      .  (8) 
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2.  3    ICP  Repaired  Inventory  Distribution  and  Order  Rules 

The  repair  facility  may  not  be  able  to  start  on  a  full  ICP  order  due 
to  an  insufficient  NRFI  inventory  of  reparable  carcasses  on  hand.      The 

point  in  time  when  a  new  repair  order,     Q  ,    will  affect  the  inventory 

R  +   1 

at  the  ordering  activity  depends  upon  the  leadtime  and  the  amount  of 

shortage  at  the  repair  facility  receiving  the  order.      As  pointed  out  by 

Hanssman   [  5]  ,    this  shortage  is  characterized  by  the  average  time    K 

it  takes  the  repair  facility  to  release  the  full  amount  of  the  order  quantity 

Q  .     Referring  to  Figure  1  ,    the  actual  output  by  the  repair  facility 

R  +   1 

into    q  will  be  equal  to  the  total  unfilled  orders  or  the  total  supply 

on  hand,    whichever  is  smaller.      It  is  assumed  that   K      is  an  integer 
multiple  of  the  unit  period.      Therefore,    the  first    I      inventory  level  that 
the  ICP  can  control  is 

I^  (R    +    K^     +     1)      . 

Now, 


I^  (R    +    K^    +    1)    =    1^(0)    +    O^    +    Q^    +    ■■'    +    Q^   +    Q^^  ^ 


(^^1    '    ■'•     +    ^^R+Kl    +    1    ^    -^'^ 


Let   H      be  defined  as  the  high  limit  of  the  repaired  item  system  which  is 
a  fixed  value  to  be  determined  by  management.      This  limit  is  the 
maximum  level  of  items  represented  by  the  ICP  repaired  item  inventory 
and  outstanding  repair  order  quantities.      (The  repair  facility  inventory 
is  not  considered  within  this  limit.  )     Each  period  an  order  is  placed  to 
bring  the  repaired  item  system  back  up  to  this  level.      It  follows  that 
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1=1 

or  (10) 

R 
Q  =    new  order    =    H„     -     1.(0)     -      E     Q. 

X\.     +       1  C  L,  .         ,  1 

1=1 


Based  on  (10)  ,    we  can  write   (9)    as  follows: 

r+kVi 

I     (R    +    K      +     1)    =    H       -  r         rX.         .  (11) 

1=1 


From   (9)  ,     it  can  be  shown   (similar  to  the  Appendix  A   derivation) 
that    I      is  a  normally  distributed  random  variable  with  parameters 

^^      =    H^     -    (R    +    K^     +     l)r^t^  (12) 

and 


a^      =    (R    +    K^    +     1)  r^a^^      ,  (13) 


where  a  negative  inventory  is  considered  to  be  the  outstanding  backorder 
quantity  when  the  on-hand  inventory  level  is  zero. 

By  the  method  shown  in  the  equation  (7)    derivation,    the  value  of  the 
mean  parameter   jj,       can  be  written  as  a  function  of  a  desired  protection 


level  a  as  follows: 


^^    =    N(a)   (R    +    K^    +    1)^^^   ra^        .  (14) 


From  equations    (12)    and   (14)  ,    the  value  of  the  high  limit  of  the 
repair  system   H      can  be  expressed  as  a  function  of  the  desired  protection 
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level   a   and  the  parameters   ij,       and  a 


^2     -' 


"N(a)  (R  +   K     +    1)^^     a^    +    (R  +    kW    l)^^  ] 


(15) 


2.4    ICP  Combined  Inventory  Distribution 

Having  determined  the  distributions  of   I      and    I     ,     the  parameters 
of  the  combined  inventory  distribution,     1=1     +    I     ,     can  now  be 
determined.      Due  to  the  common  demand  distribution  from  which    I 
and   I      were  derived,    a  dependent  relationship  exists  between    I      and    I     . 
Under  these  conditions,    the  parameters  of  the  combined  inventory 
distribution,    I,    are  defined   [6]    as  follows: 

and 

a^    =    0^    +    o^      +    2COV(I^,    I^)      .  (17) 

Under  the  assumption  that  the  joint  distribution  is  normal,    it 
follows  that  the  combined  inventory  distribution,    I,     is  also  normally 
distributed  with  the  parameters  as  shown  above.      An  alternative  form 
of  the  mean  parameter  can  be  obtained  from  (7)  and  (14)    and  the  as- 
sumption that  the  same  protection  level   a   will  be  set  for  both  the    I 
and   I      inventories. 

M,     =     N(a)a-^   (  (L+   1)^^^  (1    -  r)    +     (R   +    K^    +    l)^^^r)     .        (18) 

From   (5)  ,    ( 1  3)  ,    and  the  derivation  of   GOV  (I,    I    )    in  Appendix  B, 
the  variance  parameter  can  be  expressed  as  follows: 
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a^     =    (  (L  +   1)  (1    -  r)^    +     (R   +    K^    +    l)r^ 

+     2r  (1   -  r)  (R  +    K^    +    1)  )o^        .  (19) 

The  above  variance  equation  only  holds  if  it  is  assumed  that  leadtime 
relationships  satisfy 

L    ^    R    +    K 
Further,    although  this  equation  gives  insight  as  to  the  effect  of  the 
various  system  parameters  upon  the  combined  inventory  variance,    it 
cannot  be  used  to  develop  a  single  protection  level  relationship  unless 
the  separate    I      and   I      inventory  assumption  is  relaxed.      To  relate  a 
protection  level  to  this  term    (vice  separately  for    I      and    I    )   would 
imply  that  material  from  one  inventory  could  satisfy  demand  being 
experienced  by  the  other  inventory.     In  terms  of  practical  application 
of  the  model,   this  interchange  would  be  a  highly  desirable  extension; 
it  will  be  discussed  in  subsection  2.  5  . 

2.  5   Repair  Facility  Inventory  Distribution 

Similar  to  the    I     demand  derivation  in  Appendix  A,    it  can  be  shown 

that  the  random  demand,    rX  ,    charged  against  the    I      inventory  in  a 

2        2 
given  period  is  normally  distributed  with  parameters    (rjo,      ,    r     j       )  . 

Based  on  the  repair  order  rules  developed  in  subsection  2.  2,   the 

corresponding  ICP  repair  order  determined  at  the  end  of  each  period  is 

equivalent  to  the  repaired  item  demand  during  that  period.     It  follows 

that  the  ICP  repair  orders  can  be  thought  of  as  being  generated  from 

the    I      demand  distribution  given  above.      Further,    the  input  into  the 
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repair  facility  each  period   (number  of  reparable  carcasses  returned) 
was  assumed  to  be  generated  from  the  same  demand  distribution 
K  periods  earlier.     Assuming  that   K   is  greater  than  one  period,    the 
repair  facility's  net  inventory  position  each  period  is  a  random  vari- 
able determined  by  the  difference  between  these  two  independent 
(K  periods  apart)    normal  random  variables.      Given  these  conditions, 
it  follows  that  the  repair  facility's  inventory    V    is  normally  distributed 
with  parameters  defined  as  follows: 

M.^     =     0  (20) 

and 

2  2        2 

where  the  net  value  of  the  inventory  takes  on  negative  values    (backlog) 
when  the  ICP  repair  order  exceeds  the  on-hand  physical  inventory. 

In  subsequent  sections  where  inventory  carrying  costs  are  considered, 
it  becomes  necessary  to  determiiie  the  expected  overage  so  that  the 
average  holding  cost  of  the  on-hand  physical  inventory  can  be  estimated. 
This  value  is  found  in  the  following  manner. 

Defining  the  overage  as 

V  if      V   s    0 

overage    =    i  , 

0  otherwise 

the  expected  overage  is 
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E  [overage]     =   ——i^ V   exp  f    -  —    j 


Y^V   a 


2 

^  ^    dV 


(22) 


^V  "^X 


The  expected  shortage  value  can  be  determined  in  a  similar  manner. 
The  average  number  of  periods,    K    ,    it  takes  the  repair  facility  to 
release  the  full  amount  of  a  repair  order  can  be  determined  from  the 
expected  shortage  value. 

Defining  the  shortage  as 


shortage    =    J 


V         if      V    ^    0 


0  otherwise 


the  expected  shortage  is 

0 
E   [shortage]     =     -      -— J        Vexp^-  —^ 


^^"v 

J    ' 

-  00 

''v 

"'x 

^% 


(23) 


V2^ 


Since    r  |j,       is  the  mean  input  into  the  repair  facility  each  period,   the 
expected  number  of  periods  required  to  accumulate  sufficient  material 
to  fill  the  expected  shortage  can  be  represented  approximately  by 


E  [ shortage]  o 


""^X  '^^X 
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Obviously,    the  preceding  quotient  will  include  fractions  of  periods. 
However,    in  keeping  with  the  subsection  2,  3    derivation,    it  becomes 
necessary  to  treat  a  fraction  of  a  period  as  an  additional  whole  period. 
This  represents  a  conservative  approach  in  opposition  to  the  alternative 
of  dropping  the  fraction,    i.  e.  ,     K      will  be  a  fraction  of  a  period  larger 
than  the  expected  delay  time  vice  a  fraction  smaller.      Accordingly,     K 
will  be  defined  as  follows: 


1          r             X  ~1 

K       =    I    —= +     1       .  (24) 


where    [    ]    is  defined  as  the  greatest  integer  in  the  enclosed  value. 

2.  6    Operating  Cost  Considerations 

As  pointed  out  in  [  2]  ,    there  are  numerous  costs  associated  with 
the  management  of  a  reparable  inventory  system.      Since  the  construction 
of  a  mathematical  model  of  an  inventory  system  is  motivated  by  a  desire 
to  improve  the  operating  rules  for  controlling  a  particular  inventory 
system,    the  pertinent  costs  are  those  which  are  influenced  by  the 
operating  doctrine.      It  follows  that  costs  that  are  independent  of  the 
operating  doctrine  need  not  be  included  in  the  analysis.      The  periodic 
review  costs  considered  herein  have  been  so  restricted  and,    wherever 
possible,    have  been  grouped  under  one  symbol  to  reduce  the  model's 
notation  complexity.      The  costs  are  defined  as  follows: 

(a)      Purchase  Order  Cost   (A     )  .      The  fixed  costs  associated  with 
(1)    reviewing  the  purchase  order  records  and   (2)    preparing  a 
purchase  order. 

23 


(b)  Repair  Order  Cost   (A     )  .      The  fixed  costs  associated  with 

R 

(1)    reviewing  the  repair  order  records,    (2)    preparing  a  repair 
order,    and   (3)    initiating  batch  repair  action  at  the  repair 
facility  in  response  to  a  repair  order   (regardless  of  size)  . 

(c)  ICP  Holding  Cost   (h)  .      The  cost  per  unit  per  basic  review 
period  to  hold  a  RFI  item  in  the  ICP  inventory.      This  cost  is 
normally  a  function  of  item  purchase  cost  and  an  estimated 
holding  rate. 

(d)  Repair  Holding  Cost   (h     )  .      The  cost  per  unit  per  basic  review 
period  to  hold  a  NRFI  item  at  the  repair  facility.      This  cost  is 
normally  a  function  of  the  NRFI  item  value  and  an  estimated 
holding  rate. 

Since  the  net  inventory  distributions  are  subsequently  used  to 
determine  the  value  of  the  holding  costs,    it  is  emphasized  that  these 
distributions  were  derived  by  treating  the  random  demand  each  period 
in  a  discrete  manner;    i.  e.  ,    as  if  the  total  period  demand  were  received 
at  a  point  in  time  within  the  period.      Consequently,    a  holding  cost  is 
incurred  only  if  the  period's  on-hand  inventory  plus  receipts  exceed 
period  demand.      The  period  holding  cost  is  then  proportional  to  the 
excess  stock  on  hand. 

A  specific  shortage  cost  has  not  been  included  in  the  analysis; 
however,    a  shortage  cost  is  implied  by  the  decision-maker's  choice  of 
protection  level.      In  this  regard,    it  is  envisioned  that  management  would 
be  presented  with  yearly  cost  estimates  for  various  protection  levels  and 
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that  the  final  protection  level  decision  would  be  influenced  by  the 
associated  costs.     Accordingly,    the  primary  objective  of  deriving  a  cost 
expression  for  the  basic  model  is  to  provide  a  means  of  relating  dif- 
ferent protection  levels  to  the  expected  annual  operating  costs  corre- 
sponding to  these  levels  of  protection. 

The  analysis  will  be  based  on  the  periodic  structure  of  the  basic 
model  presented  in  subsection  2.  1,    where  the  purchase  and  repair 
orders  are  initiated  at  the  end  of  each  period.      Defining    T    as  the  basic 

period  length  in  year  units,    the  average  annual  purchase  order  cost 

A^  A^ 

is   — — —   .       Similarly,    the  average  annual  repair  order  cost  is 


—  J ,  _ .         jj _  _  —^  __ _      .,  — 

The  expected  period  cost  of  holding  inventory  at  the  ICP  is  found 
by  multiplying  the  ICP  holding  cost   (h)    times  the  expected  value  of  the 
net  inventory  at  the  ICP.     Assuming  that  management  will  set  the 
protection  level  for  the    I      and   I      net  inventories  sufficiently  high  so 
that  backorders  are  incurred  only  in  small  quantities,    the  expected 
values  of  the  net  inventories  will  very  closely  approximate  the  expected 
values  of  the  on-hand  physical  inventories.     Accordingly,    the  mean  value   |j, 
(where    |j,    =   p, ,    +   (j,^  )    of  the  combined  ICP  inventory,    I,     will  be  used  to 
approximate  the  expected  value  of  the  physical  inventory.      The  expected 

holding  cost  per  period  is,    therefore,     hjj,     and  the  expected  yearly  cost 

hn, 
IS     — '^ 


T 

The  expected  period  cost  of  holding  inventory  at  the  repair  facility 

is  found  by  multiplying  the  repair  holding  cost   (h     )    times  the  expected 

value  of  the  physical  inventory  at  the  repair  facility.     Here  it  cannot  be 
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assumed  that  backorders  are  incurred  in  small  quantities.     As  a  result 
of  the  derivation  in  subsection  2,  5,    a  backorder  or  overage  position  is 
equally  likely  at  the  repair  facility.     Using  the  expected  value  of  the 
physical  inventory   (expected  overage)    developed  in  subsection  2.  5,    the 
expected  period  cost  of  holding  inventory  at  the  repair  facility  is 


V^ 


The  expected  annual  cost  of  holding  inventory  at  the  repair  facility  is 
then 


All  the  terms  to  be  considered  in  the  expected  annual  cost  expression 
of  the  basic  model  have  now  been  evaluated. 


=   tC^p^^r^i^^  ^  T7=^)    •  "'> 


Annual  Cost 

Or,    replacing  |j,    with  its  equivalent   (18)    form,    the  annual  cost  can  be 

expressed  as  follows: 

1 

1  y 

Annual  C 


ost    =    -1-    (^Ap  +   A^+   hN(a)  a^  ((L  +    1)^   (1    -   r) 

1               , 
+    (R   +    K      +1)      r)    +       ^ ^  ^ 


1  ,  (27) 


1         ,>2     .  V     "X 

Equation  (27)  relates  the  annual  operating  costs  with  the  protection 
level   a.     As    a   is  increased,    the  value  of  N(a)   will  increase,    causing 
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the  annual  operating  cost  to  increase.     Since  the  basic  model  assumed 
a  given  review  period  length,    the  operating  cost  can  only  be  varied  by 
varying  the  protection  level.      Therefore,    once  the  basic  review  period 
has  been  established,    the  decision-naaker 's  final  protection  level  choice 
will  be  governed  by  the  associated  annual  operating  costs  he  is  willing 
to  accept  or,    in  the  case  of  a  budget  constraint,    the  annual  operating 
funds  available. 
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3.       EXTENSION  OF  THE  BASIC  MODEL 

The  structure  of  the  basic  model  was  based  upon  a  given  review 
period.     Review  and  order  action  took  place  at  the  end  of  this  unit 
period  and  ordered  quantities  were  received  at  the  beginning  of  the 
period.      The  following  extension  of  the  model  to  evaluate  the  effect  of 
extending  the  ICP  repaired-item  review  period  and/or  the  purchased- 
item  inventory  review  period  will  maintain  the  same  relative  time  point 
relationships  between  these  events  within  the  extended  periods.     This 
requires  that  only  those  time  periods  that  are  integer  multiples  of  the 
basic  period  and  are  evenly  divisible  into  the  original  number  of  lead- 
time  periods  can  be  considered  in  the  analysis. 

Let    T      be  defined  as  an  extended  time  period  with  length  equal  to 
an  integer  multiple  of  the  basic  time  period  such  that     L   s    T      ^1 

and    -— —     is  an  integer.      Let    T      be  similarly  defined  as  an  integer 

'  R 

multiple  of  the  basic  time  period  such  that    R   s:    T      s    1     and    -= —     is 

an  integer.      The  denaand  over  these  extended  periods  is  the  sum  of  the 

# 
independent  basic  period  demands  contained  within  the  extended  interval. 

It  follows  that  the  demand  over  the  extended  periods  remains  normally 

distributed.      The  mean  and  variance  parameters  of  the  extended  period 

demand  distributions  can  be  obtained  by  scaling  the  corresponding  basic 

nnodel's    I      and   I      demand  parameters    (mean  and  variance)   by  the 

integer  values  of   T      or    T     .      The  length  of  the  purchase  and  repair 

order  leadtimes  in  terms  of  the  extended  periods  are  represented  by 


_L  R 

T 


and    -7=r—    ,    respectively. 


1  2 
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Using  the  preceding  definitions,    all  the  basic  model  equations  in 
section  2    can  be  rederived  in  the  sanne  manner  as  before  but  in  terms  of 
the  extended  periods    T      and    T      .     A  similar  rederivation  of  the  repair 
facility  inventory  distribution  implies  another  restriction  on  the  values 
that  can  be  taken  on  by    T     .      Recalling  that  the  subsection  2.  5    derivation 
assumed  that  the  repair  facility  input  was  independent  of  the  ICP  repaired 
item  demand  during  the  previous  period,    it  becomes  necessary  to  further 
restrict    T      to  values  less  than  the  turn-around  time    K  .      Keeping  in 
mind  the  restrictions  on    T      and    T     ,    the  principle  equations  developed 
in  section  2    can  be  rewritten  as  functions  of   T      and    T      as  follows. 

(a)     ICP  Purchased  Item  Inventory  (I    ) 

a^    =    (L    +     T^)  (1   -  r)^a^  (28) 

i 

p,^     =     N(a)    (L    +     Xp^    (1    -  r)a^  (29) 

\_ 
if^     =    (1    -  r)    (N(a)ax  (L    +     "^i  ^  ^    +    ^^    +     '^1^^'x^  ^^°^ 


(b)      ICP  Repaired  Item  Inventory   (I    ) 

o^    =    (R    +    K^  T^    +     T^)  r^  o^  (31) 


^^    =    N(a)    (R    +    K^   T^    +     T^)^   r  a^  (32) 
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H^    =    r  (N(a)a^  (R    +    K^  T^    +    T^)^ 


+    (R    +    K^   T^    +     T^)  p,^)  (33) 


(c)     Repair  Facility  Inventory   (V) 


M-y    =     0  (34) 


2  2       2 

a„     =    2T^  r     a^  (35) 


V  2        "X 


1 


E  [overage]     =    • (36) 


^^r 


"x 


(IT    T^)  M-^ 


(d)     Expected  Annual  Operating  Cost 
Annual  Cost    =     —    (^  "Y"    +    "r"    +    hN(a)a^ 


•   (^(L   +    T^)^    (1    -  r)    +     (R   +    K^    T^   +    T^  )  ^  r^ 

i 

N(^2^'^^X         X 
+      -^^ Z= ^       )  (38) 
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The  values  of  T  and  T  that  minimize  the  annual  operating  cost 
expression  (38)  for  a  given  protection  level  (a)  can  be  determined  in 
the  following  manner: 

(1)  Hold    T      constant  at  a  value  of  one.     Solve    (38)    for  all  integer 
values  of   T      permitted  by  the  model.     Select  that  value  of   T 
which  minimizes    (38)  . 

(2)  Hold   T      at  the  value  determined  above.     Solve    (38)    for  all 
integer  values  of   T      permitted  by  the  model.     Select  that  value 
of   T      which  minimizes    (38)  . 

Using  the  above  values  of   T      and   T     ,    the  corresponding  high 

limits    H      and   H      can  be  calculated  from   (30)    and   (33)  . 
X  c^ 
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4.       EXAMPLE 

The  following  numerical  exannple  is  used  to  indicate  the  nature  of 
the  solutions  given  by  the  model.      The  parameter  values  chosen  are 
as  follows: 

T        =         1/12  year 

r  =  .9 

L  =  9  months 

R  =  4  months 

K  =  5  months 

A  =  $200 

A  =  $100 

h  =  $20   per  item  -  month 

h  =  $10   per  item  -  month 

oi  =  .  95   or    .  99 

With  these  values,    the  following  results  were  obtained: 
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9  5%  Protection   {a   =    .  9  5) 


99%  Protection  (a  =    ,  99) 


N(a)    =     1.65 


.1 


K      =     1 


N(a)    =    2.326 


.1 


K       =    1 


T       =    9    (months) 
T       =     1    (month) 
H       =     74.  8 
H^    =     229.6 
Annual  Cost   =    $5,873 


T       =9    (months) 
T       =     1    (month) 
H^     =     75.9 
H^    =     236.  5 
Annual  Cost   =    $7,  580 


The  costs  per  year  of  operating  the  inventory  system  under  the 

basic  model   (T      =    1   ,      T      =    1)    were  calculated  to  be    $7,835   and 

$9,  471    for  the  95%  and  99%  protection  levels,    respectively.      The 

combination  of   T,     =9    and    T      =    1    produced  the  minimum  annual  costs 
1  2  ^ 

under  both  the  selected  protection  levels  as  shown  above.      The  other 
combinations  permitted  by  the  model,    i.  e.  ,     T      =    1  ,    3  ,    9    , 
T      =    1  ,    2  ,    4  ,     generated  higher  annual  costs. 

The  model  is  highly  sensitive  to  the  standard  deviation  of  demand, 
a^  .      For  example,    if  a„   was  actually   8    vice    4  units  as  used  in  the 
above  example,    the  99%  protection  level  would  be  reduced  to   87% 
(holding  the  cost  and  other  operating  variables  constant)  .     Retaining  the 
same  review  period  combination     T      =    9    and   T      -    1   ,     this  sensitivity 
can  be  illustrated  in  another  manner.      If  a^   were  raised  from   4   to    8 
units,    the  same  99%  protection  level  computation  would  cause  the  ex- 
pected annual  operating  cost  to  almost  double  from    $7,  580   to    $14,  891   . 


33 


The  sensitivity  of  the  model  to  the  repair  and  purchase  order 

leadtimes,    R  and   L,    depends  upon  the  values  of  o       and   r  .     As   a 

X.  X 

increases,   the  model  becomes  more  sensitive  to  the  leadtimes  which 
are  scaled  by   r    and  (1    -  r)  ,    respectively.      For  example,    if  the  repair 
order  leadtime  for  the  99%  protection  example  was  reduced  from  4  to 
2   months,   the  expected  annual  operating  cost  would  be  reduced  only 
$903   from    $7,580   to  $6,677.     However,    if  a„   were    8    vice   4   units, 
this  same  reduction  in  leadtime  would  cause  the  expected  annual 
operating  cost  to  be  reduced   $3,004   from    $14,891    to    $11,887  . 

The  selection  of  the  protection  levels    (9  5%  and  99%)    for  the  ex- 
amples was  arbitrary.     A  family  of  solutions  could  be  obtained  for 
various  protection  levels  and  presented  to  management  for  the  selection 
of  that  protection  level  that  best  met  the  funding  constraint  for  the 
particular  item  in  question. 

The  fact  that  the  example's  minimum  cost  was  found  at    T      =9 
points  up  the  limitation  of  the  extended  model.      This  value  represents 
an  end-point  solution  in  that  this  is  the  maximum  value  that  can  be  taken 
on  by    T      within  the  extended  model  structure.      The  possibility  that 
additional  savings  might  be  realized  by  extending  the  review  period 
length  even  further  remains  unresolved.     Looking  in  the  other  direction, 
it  is  possible  to  evaluate  the  effects  of  period  lengths  smaller  than  the 
example's  basic  unit  period    T   by  equating    T   to  a  smaller  period,    i.  e.  , 
a  week  vice  a  month,    and  using  correspondingly  smaller  demand  data. 
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Nevertheless,    the  relatively  few  discrete  points  at  which  annual 
operating  costs  can  be  evaluated  remains  a  limiting  characteristic  of 
the  model. 
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5.       CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  the  assumptions  and  the  evaluation  restrictions  of 
section   3  ,    the  model's  scope  of  application  is  necessarily  confined  to 
particular  cases.      Within  this  scope,   the  equations  provide  a  means  of 
obtaining  an  understanding  of  the  interactions  of  the  principle  parameters 
and  how  they  affect  the  inventory  positions  and  costs  associated  with 
operating  a  system  of  this  nature.      The  model's  equations  indicate  the 
sensitivity  of  the  reparable -item  system  inventories  to  order  leadtimes, 
return  rates,   mean  demand,    and  protection  levels.     More  significantly, 
the  equations  highlight  the  high  sensitivity  of  the  system  to  the  demand 
variance  parameter. 

Possible  courses  of  action  to  improve  the  model  of  the  periodic 
review  system  addressed  by  this  paper  include  the  following: 

(1)      Define  an  acceptable  method  and/or  criterion  whereby  material 
from  either  one  of  the  two  ICP  inventories  treated  by  this  paper  could  be 
used  to  satisfy  demand  experienced  by  the  other.      This  could  be  permitted 
in  the  present  basic  model  if  one  accepts  the  rationale  that  the  ordering 
rules  are  based  upon    "paper"    inventory  positions  and  that  the  physical 
inventories  could  be  co-mingled.      This  would  permit  the  use  of  the  basic 
nnodel's  combined  inventory  variance  term   (19)    in  an  overall  protection 
level  computation.     However,    under  the  present  structure,    the  covariance 
between  the  two  inventories  would  be  difficult  to  evaluate  in  the  extended 
model  cases  due  to  the  changing  dependence  relationships  between  repaired 
and  purchased  item  demand  over  varying  period  lengths. 
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(2)    Modify  the  periodic  review  structure  in  such  a  manner  that  the 
integer  restrictions  that  confine  the  present  model  could  be  relaxed. 
This  would  permit  finer  resolution  of  the  least  cost  review  period 
combination  and  more  exact  treatment  of  the  expected  delay  time   K 
at  the  repair  facility. 

This  paper  looked  upon  the  reparable  item  system  as  two  subsystems. 
"Optimizing"   this  type  of  structure  eventually  leads  to  combining  sub- 
optimized  parts  vice  truly  optimizing  the  whole  system.     However,    if  one 
subsystem  is  significantly  larger  than  the  other,   the  shortcomings 
normally  associated  with  suboptimization  should  be  considerably  less 
than  if  both  subsystems  are  of  equal  stature.     As  a  result  of  high  re- 
covery rates    (estimated  at   90  -95%)   being  experienced  on  reparable 
items,    the  repaired  item  subsystem  addressed  by  this  paper  would  be 
proportionally  larger  (in  terms  of  items  controlled)   than  the  purchased 
item  subsystem.     Under  these  conditions,   the  savings  that  would  be 
realized  by  optimizing  the  total  system  might  well  be  less  than  the  costs 
involved  in  developing  and  operating  under  the  more  complex  ordering 
rules  required  by  total  system  optimizing  techniques.      Clearly,    a  trade- 
off does  exist  and  a  study  to  better  establish  this  relationship  would 
serve  to  channel  future  modeling  efforts  into  the  most  promising 
direction. 
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APPENDIX  A 


The  portion  of  the  total  demand  being  satisfied  from  the  purchased- 

itemi  inventory    I      is    (1    -  r)X,    where    X   is  a  normally  distributed 

2 
random  variable  with  parameters   ^       and  o^    .     As  shown  in 

X  X 

reference    [7]  ,     the  probability  density  function  for     ( 1    -  r)  X    is  given 
by 


From   (A-  1)  ,    it  follows  that    (1    -  r)X    is  a  normially  distributed 

random  variable  with  parameters: 

mean  =    ^      (1   -  r) 

2  2 

variance    =    a^    (1    -  r) 


It  was  shown  in  subsection  2.  2  that 


L  +  1 


1 


Let 


I^     =    H^     -      r       (1    -  r)X.        .  (A -2) 

i=l 


L  +  1 


Z    =      E      (1    -  r)X. 
i.l 


Due  to  the  assumed  independence  of  the    X.'s  ,     Z    is  also  normally  dis- 
tributed with  parameters: 

IJL^.      =    (L    +     1)   P-x  (^    "  ""^ 

(A-  3) 
2  2  2 

a^      =    {I.    +     1)    a^    (1   -  r) 
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The  distribution  function    F(x)    of  the  randonn  variable    I      is 
determined  as  follows: 


F       (x)      =     P  [I       ^    x] 
1 

=     P  [H      -    Z    ^    x] 

=     P  [Z    ^    H      -    x] 

=      1    -    F^(H^    -    X) 

H      -    X 
=      1    -    J  i^{^)  dz  .  (A -4) 

In  terms  of  a  standard  normal  density  function    0  (y)  ,    equation 
(A  -  4)    can  be  written  as 

^1    -  '^  -  ^^Z 


Fj    (x)     =     1     -       \  0  (y)     dy 


Therefore, 

H.     -    X 


U   (x)    = 
S  "z  -  "Z 


^z       ^  ^z  ^ 


^  [X    -    (H,     -   ^,)] 


2 


2      ^-        '-1        ^Z 

1  ^^Z 

—        e 


VI7  a^ 
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It  follows  that    I     is  normally  distributed  with  the  following 
parameters    (replacing  |j,       and  a„    with  their  equivalent   (A  -  3)    forms) 

M,^       =    H^     -    (L    +     1)  (1    -  r)p,^        , 

C^    =    (L    +     I)  (1   -  r)^  a^ 
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APPENDIX  B 
The  covariance  of   I      and    I    ,     Cov  (I    ,    I,  )  .    is  derived  as  follows: 

Covd^.    I^)    =    E[I^  I^]     -    E[I^]    E[I^]        .  (B-1) 

where     E  [X]     -    expected  value  of  X  . 

From  equations  (3)  and   (11)  ,    we  can  write 

L  +  1 


ECI^I^]     =    e[(h^     -      Z^     (1    -r)X.) 


r+kVi 


•    ("2    -      .2,         '  ^i)   ] 
1=1 


1 


H^  H^     -    H^  (L    +     1)  (1   -  r)p,^    -    H^  (R    +    K       +     l)r^J, 


X 


L  +  1  R+K^+1 


+    (1   -  r)r    E  r     L    X.  L         X.    ]       . 

i=l  i=l 


Based  on  the  condition  that    L  s   R  +   K      , 

L  +  1  R+K^  +  1      „  R+K^  +  1    ^  _  R+kVi 


r      EX.  L      X.  l    =    Er       2        X.^I+eTSL       2X.  X.  1 

i=l  i=l  1=1  i  <   j 


r+kVi         l  +  1 

+  E  [     i:  L  x.x  ] 

=    (R  +   K^    +    1)  {a^   +  ^x^    +    (R  +   kS 

•  (R+  K^    -  1)^^    +    (R+  K^  +  1)(L-R-kSm,^ 
=    (R    +    K^    +     1)  {o^    +    (L    +     l)p,^^)      . 
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Therefore, 

E  [I^   I^]     =    H^  H^  -  H^  (L  +   1)  (1    -  r)p.^  -  H^  (R  +  K^  +   l)r^, 


X 


+    (1    -  r)r    (R+  K^  +  1)    {a^   +    (L  +  l)^^)    ,      ..         (B  -  2) 


Now, 


ELIj]    E[y     =    M,^  j,^    =    (H^    -    (L+  1)    (1    -  r)^,^) 


(H^     -    (R+  kW   l)rtx^)     .  (B-3) 


Putting  the    (B  -  2)    and  (B-3)    forms  of  E  [  I     I    ]     and    E  [I    ]    E  [I    ] 
back  into  equation   (B  -  1)    and  clearing  terms  gives 

Cov(I^,    I^)     =    r(l   -  r)    (R  +  K^  +  Da^^      .  (B-4) 
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